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We investigate the response of two three-body Coulomb systems when driven by attosecond half- 
cycle pulses: The hydrogen molecular ion and the helium atom. Using very short half-cycle pulses 
(HCPs) which effectively deliver "kicks" to the electrons, we first study how a carefully chosen 
sequence of HCPs can be used to control to which of one of the two fixed atomic centers the electron 
gets re-attached. Moving from one electron in two atomic centers to two electrons in one atomic 
center we then study the double ionization from the ground state of He by a sequence of attosecond 
time-scale HCPs, with each electron receiving effectively a "kick" from each HCP. We investigate how 
the net electric field of the sequence of HCPs affects the total and differential ionization probabilities. 

PACS numbers: 03.65.Sq,32.80.Fb 



I. INTRODUCTION 

Rapid developments in the generation of attosecond 
time scale laser pulses [l| makes investigations of the in- 
teraction of multi-electron atoms with ultra-short laser 
pulses timely. The control and manipulation of Rydberg 
atoms using one or more half-cycle pulses (HCPs), with 
each HCP of duration r << r„ where T„ is the classi- 
cal electron orbital time, has been the subject of intense 
research lately [2j] . One of the tools for controlling such 
states are carefully tailored sequences of HCPs 0, y| . An 
interesting aspect of the HCPs is that atoms respond to 
them very differently than to laser pulses. The response 
of a single-electron atom to a sequence of HCPs has been 
thoroughly explored in Ref. [J]- In addition, for single- 
electron atoms, the effect that the net field of a sequence 
of HCPs has on the atomic dynamics [a, m, 0] as well 
as the remarkable control that can be achieved using a 
chirped train of HCPs @] are among the problems that 
have been studied recently. 

In the current work, we first address the effect of a 
sequence of HCPs on the electron in the H2 diatomic 
ion. The purpose of this section is to refresh concepts of 
the "kicked" one electron dynamics in one atomic center 
when driven by a sequence of HCPs 0, 0] and to explore 
new effects. We begin with the simplest scenario in which 
the two atomic centers are so far apart that the electron 
is initially localized on one atomic center. The large dis- 
tance between the fixed atomic centers allows us to make 
a connection to the response of a single electron atom to 
HCPs which has been studied intensively. In particular, 
for HCPs of picosecond duration there is a large number 
of studies for single ionization of Rydberg atoms (see for 
example [J, [a, M, [ifl ilM)- Very recently, protocols of 
generating attosecond HCP trains using strong two-color 
infrared laser pulses have been analyzed |12l | . 

The other three-body Coulomb system we explore is 
the He atom where we study the effect of an idealized 
sequence of HCPs on double ionization both when the net 
electric field vanishes (even number of kicks) and when 



it does not (odd number of kicks). 

The interaction of multi-electron systems with very 
strong and ultrashort laser pulses is still a wide-open 
problem. We perform our calculations using the Classi- 
cal Trajectory Monte Carlo (CTMC) method [ISj. Quite 
a few studies show that classical methods can be suc- 
cessful in describing the single ionization of atomic sys- 
tems when driven by ultrashort and strong laser fields, 
in very good ag reement with quantum mechanical re- 
sults [ij, \lM, 13 ■ However, open issues remain for sin- 
gle ionization of atomic systems such as the effect of the 
Coulomb interaction on the ionization process depending 
on different intensities of the field [ij, [l^, |T0| . Another 
issue concerns the range of parameters over which the 
classical techniques give the best agreement, for integral 
as well as differential probabilities, with quantum me- 
chanical studies [20, [21I, [23]. One of the goals of the 
current classical study is to initiate exploration of the 
above issues in the case of the double ionization of the 
attosecond time-scale driven He. 

The parameters we use are chosen to complement re- 
cent studies on ionization by strong ultrashort pulses 
(23I [24| . These authors have found that when the elec- 
tron is driven by two HCPs of alternating sign there is 
a very high probability for the electron to recombine to 
the ground state of H. They have derived analytic expres- 
sions, in the framework of the Magnus approximation, for 
the recombination rates for a single-electron atom (23| 
and only very recently for two electron atoms 26]. 

The current paper is structured as follows: In Section 
H we study the effect of one or two HCPs on the electron 
dynamics in the H2 diatomic ion when the two nuclear 
centers are far apart and the electron is initially localized 
in one of them. In Section HI, we obtain the total and 
double differential probabilities for the double ionization 
of the driven He atom. We explore the effect of a sequence 
of HCPs when the net electric field vanishes and when it 
does not and investigate the infiuence of the electron- 
electron interaction in the ionization process. 



II. SINGLE-ELECTRON DIATOMIC 

MOLECULAR ION DRIVEN BY STRONG 

ATTOSECOND PULSES 

In recent years quite a few studies have addressed ion- 
ization of single electron atoms when driven by HCPs, 
with each HCP with duration r much smaller than the 
time the electron needs to orbit around the nucleus (Ke- 
pler period), T„. Let us describe the sequence of HCPs 
by a field of the form E(t) = YZZi {-^VEofn{t - [n - 
l)A))z, where the shape of each HCP is given by the 
function /„(t) with (n — 1)t < t < nr, and A being the 
time delay between the HCPs numbered n and n-l-1. In 
these previous studies [J, Q it has been well established 
that the effect of an A^ = 1 HCP with t < < r„ is to de- 
liver a momentum "kick" , with the energy distribution of 
the ionized electron centered around q'^/2 — \Ei\ , when 
q^ /2 — \Ei\ > 0,with Ei being the binding energy of the 
electron. If /„(t) = sin{-K{t — A)/t), which is the shape 
of the HCP we use for all subsequent calculations, then 
the strength E^) of the field is expressed in terms of the 
momentum "kick" as Eq — g7r/(2r), and the sequence 
of HCPs we use effectively delivers "kicks" of alternating 
direction. In recent studies [23|j l2J] it was shown that 
while a single {N ~ 1) HCP strips the electron from 
the hydrogen atom, a subsequent application of a second 
HCP of opposite sign immediately after the end of the 
first one results in re-attachment of the ionized electron 
with probability close to one. 

Studying the above concepts in a new context we now 
investigate the effect of a sequence of iV = 2 HCPs on the 
single electron of the H^ ion. We assume that initially 
the two nuclei are along the z-axis at i? = 40 a.u. and 
the electron is localized on the left nucleus. Since the 
distance R of the nuclear centers is large it is a very 
good approximation to consider the initial state of the 
single electron as that of a hydrogen atom. We take 
as our initial state the 2s state of the hydrogen atom 
(left nucleus) and consider the two nuclei fixed in space. 
The initial distribution of the 2s state is given by the 
microcanonical distribution jl3| p(7): 



p{^) ^ .yV X S{Ei + h) 



(1) 



with the normalization constant ,A^ and /i — 1/8 a.u. 
namely the ionization energy of the 2s electron. 



A. Single electron ionization in the J/j w^hen 
driven by an TV = 1 HCP 

We first study the effect of an iV = 1 HCP on the 
diatomic ion using the full-three body Hamiltonian with 
nuclei fixed in space. 



H = pV2 - l/|r - Ril - l/|r - R2I - r ■ E(t), (2) 



with Ri , R2 the position vectors of the two nuclear cen- 
ters and E(i) the field specified above polarized along 
the direction of the molecular axis z. The probability 
for ionization from the diatomic ion as a function of the 
strength of the momentum transfer q is shown in Fig. [T] 
for a HCP with t = 3 a.u. Since the nuclei are so far 
apart it is of course to be expected that our results in 
Fig. [T] are exactly the results one would obtain from the 
single ionization of hydrogen from the 2s state, as is in- 
deed the case [24J. In our classical calculation with the 
initial distribution considered in Eq. ([T]) the agreement 
with quantum mechanical results is better the larger the 
strength of the field is and thus the higher it is from the 
threshold field strength corresponding to over the barrier 
ionization. Let us also emphasize that for all the results 
presented the full Hamiltonian of Eq. ^ is used for the 
propagation in time and our reference to "kicks" is only 
an interpretation of the accurate results obtained with 
all interactions accounted for. 



B. Controlling the atomic center of electron 
re-attachment in H2 using N — 2 HCPs 

Next, we consider the case of A^ = 2 HCPs, with q — 2 
a.u. and t — 3 a.u., where we vary the time delay A 
between the two HCPs. What we are effectively doing is 
first "kicking" the electron from left to right with a mo- 
mentum transfer of g = 2 a.u. and then "kicking" it in 
the opposite direction with a delay A. At the end of the 
two HCPs we compute the electrons that remain bound, 
that is, the trajectories for which the energy of the elec- 
tron given from Eq. ^ is negative. The electrons are 
bound either to the left or right nuclear center which we 
specify by checking whether z < or 2: > respectively, 
where z is the coordinate of the electron along the molec- 
ular axis, with the of the axis being at the center of the 
two nuclei. The results are shown in Fig. [21 As expected, 
when the two HCPs which have opposite direction and 
equal strength are delivered one immediately after the 
other the electron remains almost completely bound to 
the left nucleus and the probability to be bound to the 
right one is zero. However, as the delay between the 
two pulses increases the probability to find the electron 
bound to the right nucleus significantly increases. This is 
reasonable since the longer the delay time, the more time 
the ionized electron from the left nuclear center needs to 
travel further towards the right nucleus. When the sec- 
ond "kick" is received in the opposite direction it reduces 
the kinetic energy of the electron ionized from the left nu- 
cleus, and depending on the electron's position from the 
right nuclear center it can result in re-attachment of the 
electron to the right nucleus. Note that the maximum 
of the re-attachment probability takes place at a A very 
close to the time it takes for the electron with a momen- 
tum approximately q = 2 a.u. to travel from the left 
to the right nuclear center R/q = 40/2 = 20 a.u. Note 
that our three-dimensional calculation indicates that the 



probability for the electron to be re-attached when one 
uses HCPs is very large which is to be expected since the 
ionization probability from the first HCP is very large 
launching with high probability electrons to the contin- 
uum and thus resulting in higher recombination probabil- 
ity to the other nuclear center. With rapid experimental 
advances, this method could be used in the future for 
determining distances between atomic centers. For re- 
lated work on how to use attosecond pulses (not HCPs) 
to observe eflFects similar to the ones discussed above in 
diatomic ions see ref. Ha, SI, Ell ■ 



III. TWO ELECTRONS DRIVEN BY STRONG 
ATTOSECOND PULSES 

We now turn from single electron dynamics in one 
atomic center to two-electron dynamics in two atomic 
centers. In the case of a driven two-electron atom the 
Hamiltonian is 



H 



Pi 
2 



2 ri r2 



1 



|ri -ral 



(ri+r2) •£(*). (3) 



In what follows we use as an electric field the same lin- 
early polarized sequence of HCPs as for the case of the 
one electron atom in the previous section with the time 
delay being fixed now and equal to r (i.e., effectively we 
have a sine pulse). 



Initial phase space distribution for the 
two-electron atom 



Ze//2 for electron 1 and 2, respectively. We then use the 
initial conditions we have just generated to obtain the 
total energy of the two electron atom from Eq. ([3|) by 
setting the field equal to zero. The initial total energy 
distribution of the two electrons is shown in Fig.[3l where 
we see that there is a long tail on both sides of the ideal 
value of -2.9 a.u. We cut the tail off by introducing two 
parameters Emm, Emax such that [3l| : 



/|-; EpiE)dE 



E„ 



(5) 



with Ep the most probable energy and p{E) the energy 
distribution. In our calculations Emax = —2.51. ii) In 
addition, we evolve the two-electron atom freely (i.e., 
without an external field and using the full three-body 
Coulomb Hamiltonian) and discard the trajectories for 
which the total two electron energy becomes positive dur- 
ing the field-free propagation. We freely propagate the 
system for times twice the Kepler period, with the latter 
being approximately 2.2 a.u. for the He ground state. In 
this latter step we find that the fraction of trajectories 
labeled as "auto-ionizing" and thus discarded in our sim- 
ulation is small. After these steps, we obtain an initial 
state distribution which has the radial and momentum 
distribution shown in Fig. S) Finally, we also find that 
the field-free evolution of the initial momentum distri- 
bution changes little for times comparable to the Kepler 
period (see Fig. [5]) and we can thus consider our ensem- 
ble of initial conditions as approximately stable for all 
practical considerations. 



The initial phase space density p{"f) in our classical 
calculation of the double ionization of He is given by a 
product of microcanonical distributions 



P{l) = ^S{ei + Ii)6{e2 + h), 



(4) 



with normalization constant ^. In the case of He we 
account for the electron-electron repulsion in the initial 
state indirectly through the use of effective charges [30] ■ 
In the following we present results for two sets of effective 
charges: a) h ~ I2 = -Z^ef f/2 where Z^ff — 27/16 and b) 
/i = Z\r /2 and I2 — ZKr/2 with Ii and I2 being the 
ionization potentials for the two electrons in the Is^ state 
of He, i.e., Ii = 2 and I2 = 0.9 a.u. The latter choice of 
effective charges accounts better for the electron-electron 
repulsion in the initial state. 

A few more remarks on the construction of the initial 
state are in order: Our initial state is not an eigenstate 
of the driven two-electron Hamiltonian in Eq. ([3|), and 
it can thus "auto-ionize" even without the presence of 
the external field. We avoid the latter problem in the 
following way: i) We first generate the initial conditions 
for each electron independently using the microcanonical 
distribution for a one electron atom with charge •Z^e//i or 



B. Computation of the doubly ionizing trajectories 

We determine doubly ionizing trajectories due to the 
HCP's as follows: when the pulse is switched off we check 
the total energy of the three-body Coulomb system. If 
the energy is negative then the driving field has not trans- 
fered enough energy to cause double ionization and we 
thus label those trajectories as non-doubly ionizing tra- 
jectories. If the energy is greater than zero then it can be 
the case that these trajectories will lead to single or dou- 
ble ionization. To decide if these latter trajectories are 
singly or doubly ionized we continue to monitor them 
in time by propagating the freely evolving three-body 
Coulomb system until the asymptotic regime is reached. 
If at that time e; = pf/2 — Z/ri, with Z = 2, for each 
electron are both positive we label these trajectories as 
double ionizing otherwise as non double ionizing. Note, 
that the above described method does not allow us to 
separate the singly ionized from the bound trajectories. 
The reason is that with the process described above the 
ti — p1/2 — ZJTi reach constant values only asymptoti- 
cally but when the system is propagated for so long the 
classical nature of our calculation can cause bound trajec- 
tories to become "auto-ionized" and thus artificially con- 



tribute to the single ionization probability. To compute 
the double ionization probability we first find the number 
of double ionizing trajectories and normalize with respect 
to the total number of trajectories that have been prop- 
agated. The single and double differential probabilities 
presented in this paper use the asymptotic values of the 
quantities plotted. We believe that even though we ac- 
count for the electron-electron correlation in the initial 
state through effective charges this simple initial state 
captures accurately the essential final correlations as has 
been shown in other cases fSOll. 



Single and double differential probabilities for 
double ionization 



Following the procedure previously described, we con- 
struct the initial state distribution both for equal and 
different effective charges and find the trajectories that 
doubly ionize after A^ = 1 HCP. We find that both initial 
state distributions yield similar results both for the en- 
ergy as well as the inter-electronic angular distribution, 
as can be seen in Fig. [6] and Fig. [T] In the following 
we present results for the initial state distribution that 
corresponds to Ji 7^/2, In addition, for the single dif- 
ferential probabilities presented in this paper for /i 7^ I2 
we take the average of the differential probabilities of the 
two electrons unless otherwise specified. 

In Fig. [51 we present results for the energy distribu- 
tion for odd or even number of HCPs for q — i a.u. and 
T = 0.1 a.u., with r much smaller than the Kepler period 
T„ of an electron in a hydrogenic atom, T„ = 2Trn^ /Z'^^^r. 
Let us note that while an even number of HCPs corre- 
sponds to a zero time integral of the complete pulse and 
it can thus be produced, the odd number of HCPs cor- 
responds to a non-zero time integral of the pulse and 
can thus not be produced ^^. However, what can be 
produced instead of an A'' = 1 HCP is a very short and 
strong first half cycle followed by a much weaker and 
much longer second half to compensate Q. While the 
interaction of each of the electrons with the laser pulse 
is very strong, the electron-electron interaction as well as 
the interaction of each of the electrons with the nucleus 
is not negligible. If the latter were of no importance, 
double ionization would not take place, since an even 
number of HCPs transfers zero net momentum. This has 
already been noted for the ionization of the driven hy- 
drogen, see Ref. 24]. There it was pointed out that 
one can think of the process as a sequence of 6 kicks, at 
times equal with an integer multiple of t/2, with oppo- 
site signs (HCPs of opposite sign) where the interaction 
of the electron with the nucleus in between the kicks can 
not be neglected. For an even number of HCPs it is no 
surprise that on average the momentum transfer is small 
and the kinetic energy acquired is less than the bind- 
ing energy of the two-electron atom, Ei, with the energy 
distribution peaking at energies close to zero. However, 
when an odd number of HCPs are applied then the net 



momentum transfer is q and since the kinetic energy for 
each electron is greater than the electron's binding en- 
ergy, g^/2 > \Ei\/2, one expects the energy distribution 
to peak around g^/2 — |£'i|/2 w 3.1 a.u. if the electron- 
electron interaction plays no role. From Fig. [8] we see 
that the electron distribution peaks at an energy smaller 
than 3.1 a.u. The latter effect can not be due to our 
approximate initial state since this latter approximation 
can only cause the energy distribution to peak at a higher 
energy (effectively, smaller Ei). As we increase the num- 
ber of odd HCPs the shift of the energy distribution to 
smaller energies becomes more substantial, although one 
can clearly see a "shoulder" structure around i;^ = 3.1 
a.u. This shift to smaller energies for an increased num- 
ber of odd HCPs must be due to the increasing signifi- 
cance of the electron-electron repulsion in the dynamics 
of the doubly ionizing trajectories for the application of 
longer pulses. Indeed, for an even number of HCPs where 
the electron-electron interaction is more important (see 
also next paragraph) the energy distribution is peaked 
around small energies. 



In Fig. [HI we show the distribution of the inter- 
electronic angles for odd as well as even numbers of 
HCPs. We see that for an odd number the distribution 
peaks at approximately 45° while for an even number 
of HCPs its peak shifts from larger to smaller angles as 
the number of HCPs increases. The fact that the inter- 
electronic angle between electrons that have received an 
even number of kicks is larger than the angle between 
those which have been kicked an odd number of times is 
to be expected: In the former case the electrons are much 
slower (compare the right and left panels in Fig. [8]) and 
thus the electron-electron repulsion is much more pro- 
nounced as compared to the electrons that escape after 
an odd number of cycles. This will be further illustrated 
in the following when we present single and double dif- 
ferential probabilities in angle. 



In Table I we present the probability for double ion- 
ization as a function of the number of HCPs. We find 
that the double ionization probability follows the same 
pattern as the single ionization of H [2J| . Namely, as the 
number of even HCPs increases the probability for dou- 
ble ionization increases whereas the probability for dou- 
ble ionization decreases with increasing number of odd 
HCPs (Table I shows that the probability for both elec- 
trons to remain bound after a large number of even cy- 
cles is significant.) So as for the case of a single electron 
atom for the case of two electron atoms the electrons are 
"stripped" from the atom when driven by an odd number 
of HCPs while an even number of HCPs causes the two 
electrons to re-attach to the atomic center [25|. Also, in 
Fig. [in] we show how the double ionization probability 
changes as a function of the momentum transfer q for 
the case of TV = 8 and iV = 9 HCPS. 
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In addition, we show in Fig. \TT\ and Fig. [12] the prob- 
ability density of the momentum along the z-axis (see 
Ref. 33] for the definition of the classical probability den- 
sity). It is clear that the effect of the field on each of 
the electrons is a kick received during each of the HCPs. 
We emphasize that all our results are obtained by prop- 
agating with the full three-body Hamiltonian under the 
influence of an odd or even number of HCPs of total 
duration Nt, see Eq. ([3]). Referring to "kicks" is only 
an interpretation consistent with the results obtained ac- 
counting for all Coulomb interactions and accounting for 
the effect of the finite total duration of the external field 
applied. The probability density of the momentum along 
the z axis clearly illustrates that the HCPs applied are 
very strong. 

In Fig. [Island Fig.[T4|we show the double angle differ- 
ential probability densities. The angles in question refer 
to the angle that each of the electrons makes with the 
polarization axis z. We find that when the two electrons 
are driven by an odd number of HCPs they escape to the 
continuum almost antiparallel to the electric field (which 
is the direction of the force the field exerts on each of the 
electrons) with a small inter-electronic angle. When the 
number of cycles is even the two electrons escape with a 
much larger inter-electronic angle with one electron at a 
small angle with respect to the polarization axis. 

Finally, in Fig. [T5| and Fig. [T6| we plot the double dif- 
ferential probability with respect to the magnitude of the 
momenta of electrons 1 and 2. Both for the case of odd 
and even number of cycles we find that as we increase 
the number of HCPs more electrons escape with differ- 
ing momenta, that is, the probability for one fast and one 
slow electron increases with the number of HCPs. Also a 
comparison between the N = 9 and TV = 10 HCPs case 
clearly shows that the escape of the two electrons with 
different momenta is much more pronounced for the case 
of iV = 10 HCPs. In Fig. [TTj we plot the average mo- 



mentum, radius, and relative direction of the momentum 
and position, r • p/{rp), for each electron for the case of 
N = 9 and TV = 10 HCPs for all doubly ionizing tra- 
jectories. We find that at the end of the iV = 10 HCPs 
the electrons are closer to the nucleus (smaller radii) com- 
pared to the iV = 9 case. Another difference is that while 
for the case of A^ = 9 HCPs the two electrons follow the 
external field during the time that it is switched on, for 
the case of A^ = 10 HCPs the electrons strongly interact 
with the nucleus (very small value of the radii) before 
they follow the driving field. 

In this section we have presented classical results re- 
garding the double ionization of the ground state of the 
He atom with HCPs with the duration r of each HCP be- 
ing much shorter than the Kepler period. In addition, all 
our calculations for single and double differential proba- 
bilities have been performed for values of the peak of the 
external field Eq that are much higher than the threshold 
value for over the barrier ionization. So our choice of peak 
strength and duration of the pulse is such that the classi- 
cal calculations would be most reliable. We expect that 
our classical calculations accurately capture the essential 
features of the double ionizing trajectories particularly 
of those where the electrons escape with larger energies. 
Our classical results are presented for the n = 1 state 
of the He atom using parameters of the external field 
that correspond to a peak intensity of 7.8 x 10^° W/cm^ 
and an effective frequency, tt/t, of 31 a.u. which are 
well beyond current experimental capabilities. However, 
the results presented (with an appropriate scaling of the 
transfer q and the time r) should also hold for two elec- 
trons much less bound in the initial state as is the case 
for the doubly excited resonant states of He which can 
have lifetimes of picoseconds and are thus within current 
experimental capabilities. 

In summary, we have used a classical calculation to 
show how a sequence of HCPS can be used to "strip" 
and re-attach one electron to two atomic centers or two 
electrons to one atomic center. We have presented re- 
sults for the single and double differential probabilities 
for the double ionization of helium from the ground state 
when driven by HCPs of attosecond time scales. We have 
shown that for an even number of attosecond HCPs the 
effect of the nucleus on double ionization is pronounced 
and can not be neglected. We anticipate that as quan- 
tum calculations become available a direct comparison 
with our classical results will better illustrate the regime 
of validity of the classical calculations for single and dou- 
ble differential probabilities. 
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FIG. 3: Total energy distribution of the two electrons in the 
initial state without a field. 



FIG. 1; Probability for ionization of the electron from a 2s 
state of a hydrogen atom driven by an A'^ = 1 HCP with t — 3 
a.u. 
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FIG. 2: Probability for attachment in the right nuclear center 
when the electron is initially localized in the left nuclear center 
at the 2s state and is then driven by iV = 2 HCPs with g = 2 
a.u. and r = 3 a.u. 
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FIG. 4: Left figure: Radial distribution as a function of the 
radial electronic coordinate; right figure: momentum distri- 
bution as a function of the electronic momentum. Since we 
use different effective charges the two electrons are not equiv- 
alent and so we plot as a function of the radial coordinate 
of electron 1 or 2 (ri||r2) and similarly for the momentum 
distribution. 




FIG. 5: The field-free evolution in time of the initial momen- 
tum distribution. At time zero it is identical to the momen- 
tum distribution shown in Fig. |4] 
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FIG. 6: Distribution of the inter-electronic angle for a laser 
pulse with g = 3 a.u. and r — 0.1 a.u. . Results for A^ = 1 
HCP are shown: The ■ refer to Ji — I2 while the o refer to 
/i / I2. 



a, 0.15 



w 

-O 0.05 



3.1 5 10 

E ia.u.l 



FIG. 7: Energy distribution normalized over the double ion- 
ization probability for a laser pulse with q — 3 a.u. and 
r = 0.1 a.u. . Results for A = 1 HCP are shown: The ■ 
refers to Ii — I2 while the o refers to /i 7^/2. 
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FIG. 8: Energy distribution for a laser pulse with g = 3 a.u. 
and T = 0.1 a.u. Left panel: results for an odd number of 
HCPs with N = 1 (.), iV = 3 (■), iV = 5 (*), N = 7 (+), 
A'^ = 9 {o); Right panel: results for an even number of HCPs 
with N ^2 (•), N ^4 (■), iV = 6 (*), iV = 8 (+), N ^10 
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FIG. 9: Distribution of the inter-electronic angle for a laser 
pulse with q — 3 a.u. and r = 0.1 a.u. Left panel: results for 
an odd number of HCPs with A^ = 1 (•), iV = 3 (■), iV = 5 
(*), N = 7 (+), A'' = 9 (o); Right panel: results for an even 
number of HCPs with N ^ 2 {•), N = 4 (M), N ^ 6 {*), 
iV = 8 (+), iV = 10 (o). 
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FIG. 10: Double ionization probability as a function of the 
momentum transfer q for N = 8 o and A'^ = 9 o HCPs. 
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FIG. 11: Probability density of momentum along the z-axis 
for g = 3 a.u. and r = 0.1 a.u. Results for an odd number of 
HCPs with TV = 1, 3, 5, 7 and Af = 9 from left to right. 
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FIG. 12: Probability density of momentum along the z-axis 
for q = 3 a.u. and r = 0.1 a.u. Results for an even number 
of HCPs with Af = 2, 4, 6, 8 and iV = 10 from left to right. 
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FIG. 13: Double angle differential probability density for a 
laser pulse with g = 3 a.u. and r — 0.1 a.u. The angles 
Oi and 82 are the angles of the two electrons with respect to 
the direction of the field. 9i = corresponds to a direction 
parallel to the field. Results for an odd number of HCPs with 
N = 1,3,5,7 and iV = 9 from left to right. Note that the 
results have been symmetrized with respect to the diagonal 
di —62- 




FIG. 14: Double angle differential probability density for a 
laser pulse with g = 3 a.u. and t — 0.1 a.u. The angles 
611 and 612 are the angles of the two electrons with respect to 
the direction of the field. 9i = corresponds to a direction 
parallel to the field. Results or an even number of HCPs with 
iV = 2,4, 6, 8 and iV = 10 from left to right. Note that the 
results have been symmetrized with respect to the diagonal 
9i =62- (The graph for A'^ = 2 is not smooth due to the lower 
statistics) . 
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FIG. 15: Double momentum difFerential probability density 
for the momentum pi,2- The left figure corresponds to A^ = 1 
HOP and the right one corresponds to A'' = 9 HCPs. Note 
that the results have been symmetrized with respect to the 
diagonal pi — p2- 




FIG. 16: Double momentum differential probability density 
for the momentum pi,2. The left figure corresponds to A^ = 2 
HCPs, the right one corresponds to A'^ = 10 HCPs. Note 
that the results have been symmetrized with respect to the 
diagonal pi = P2- 
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FIG. 17: Left panel for A^ = 9 HCPs and right panel for 
A' = 10 HCPs: Average over all double ionizing trajectories 
of the radius of electrons 1 (black dashed line) and 2 (gray 
dashed line), of the relative direction between the position and 
the momentum of electron 1 (black solid line) and electron 2 
(grey solid line), the momenta parallel to the field for electron 
1 (black dotted line) and electron 2 (grey dotted line). The 
thin solid line represents siniji/rt). 



